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Abstract 



We investigate effects of the vector interaction on chiral and color superconducting 
(CSC) phase transitions at finite density and temperature in a simple Nambu-Jona- 
Lasinio model. It is shown that the repulsive density-density interaction coming from 
the vector term, which is present in the effective chiral models but has been omitted, 
enhances the competition between the chiral symmerty breaking (xSB) and CSC phase 
transition, and thereby makes the thermodynamic potential have a shallow minimum 
over a wide range of values of the correlated chiral and CSC order parameters. We find 
that when the vector coupling is increased, the first order transition between the xSB 
and CSC phases becomes weaker, and the coexisting phase in which both the chiral and 
color-gauge symmetry are dynamically broken comes to exisit over a wider range of the 
density and temperature. We also show that there can exist two endpoints, which are 
tricritical points in the chiral limit, along the critical line of the first order transition 
in some range of values of the vector coupling. Although our analysis is based on a 
simple model, the nontrivial interplay between the xSB and CSC phases induced by 
the vector interaction is expected to be a universal phenomenon and might give a clue 
to understanding results obtained with two-color QCD on the lattice. 
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§1. Introduction 



It is one of the central issues in hadron physics to determine the phase diagram of strongly 
interacting matter in the temperature (T)-chemical potential (/x) or T-pB plane, with pB 
being the baryonic density. In extremely hot and dense matter, the non-Abelian nature of 
QCD ensures that the colored quarks and gluons are not confined, and chiral symmetry is 
restored. Lattice simulations of QCD^^'^^ show that the QCD vacuum undergoes a chiral and 
deconfinement transition at a temperature Tc around 150—175 MeV at vanishing chemical 
potential, with the order and critical temperature being dependent on the number of active 
flavors. Although there have been several promising attempts^^"^^ to make simulations of 
lattice QCD with finite p possible, they have still not progressed enough to predict anything 
definite about the phase transition at finite pb or p. It is widely believed on the basis of 
effective theories^^^"^^) and chiral random matrix theory^^) that the chiral phase transition 
from the chiral-symmetry broken to the restored phase is first order at vanishing temperature. 
Furthermore, people believes now that the critical line of the first order chiral transition 
continues for smaller chemical potentials in the T-p plane and ends at some point with 
T — Te and p — p^, which is called the endpoint. We notice that the first-order chiral 
transition is accompanied by a jump in the baryon density. 

Recent renewed interest in color superconductivity (CS)^^^"^^^ has stimulated intensive 
studies of the QCD phase structure at finite density in the low temperature region, which in 
turn are revealing a rich phase structure of high density hadron/quark matter with CS.^^^"^^^ 
Possible relevance of CS to characteristic phenomena observed for neutron stars are being 
actively discussed. Some recent studies have also suggested that experiments on the 
Earth using heavy-ion collisions with large baryon stopping can elucidate something about 
CS in dense matter.^^^'^^^ 

The purpose of the present paper is to reveal new characteristics of the chiral to color 
superconducting (CSC) transition based on a simple effective model incorporating the vector 
interaction by focusing on the implication of the density jump accompanied by the chiral 
transition. 

Low-energy effective models^*^^"^^^'^^^'^^) are useful to study not only the chiral transition 
but also CS in dense hadronic matter. For example, chiral models of the Nambu-Jona- 
Lasinio type,^^) which can be considered simplified versions of those with an instanton- 
induced interaction, accurately describe the gross features of the T dependence of the chiral 
quark condensates of the lightest three quarks as given by lattice QCD, and predict that the 
chiral transition for 7^ is rather strongly first order at low temperatures when the vector 
interaction is absent or small.^^^'^^^'^^^ Chiral effective theories show that the gap A of CS 



3 



may become as large as 100 MeV in relatively low densities, where a phase change from the 
chiral symmetry breaking (xSB) phase to the CSC phase may also occur. ^^^'^^^ 

However, although many works on CS have been carried out with the use of effective 
models, the vector interaction, •^^^""'^^ 

Ly = -Gy{i,^^i,f, (1.1) 

has been scarcely taken into account, with the exception of very recent works. ^''^'^^^ Our 
point is that such a vector interaction is chiral invariant and naturally appears in the effective 
models derived from microscopic theories and, as we shall show, indeed can have strong effects 
on the chiral-to-CSC transition and the properties of the CSC phase. 

Although it may not be a common knowledge in the physics community, the importance 
of the vector coupling for the chiral transition is known; i.e., the vector coupling weakens 
the phase transition and moves the chiral restoration to a larger vahic of ^.36), 42), 46) 'jj-^jg 
can be intuitively understood as follows. According to thermodynamics, when two phases 
I and 11 arc in an equilibrium state, their temperatures Tj^n, pressures Pj^n and the chemical 
potentials /xiji are the same: 

Tl^Tii, Pi = Pii, = (1.2) 

If I and II are the chirally broken and restored phase with quark masses satisfying M/ > M//, 
the last equality further tells us that the chirally restored phase has a higher density than 
the broken phase, because /iiji at vanishing temperature are given by \x\ = + p'j^^ , {i = 

I, II), and hence ppj < Pfu, where pp^ is the Fermi momentum of the i-th phase. Thus it 
is seen that chiral restoration at finite density is necessarily accompanied by a density jump 
to a higher density state with a large Fermi surface, which in turn favors the formation of 
Cooper instability leading to CS. 

However, since the vector coupling includes the term {ip'-f^iljy , it gives rise to a repulsive 
energy proportional to the density squared, i.e. Gyp's f^, which is larger in the restored 
phase than in the broken phase; the vector coupling weakens and delays the phase transition 
of the chiral restoration at low temperatures. Thus one expects naturally that Cy causes 
the chiral restoration and the formation of CS to shift to higher chemical potentials, and 
may alter the nature of the transition from the xSB phase to the CSC phase drastically. 

Is it legitimate to include a vector term like (1.1) in an effective Lagrangian? First 
of all, one should notice that the instanton-anti-instanton molecule modeV^^'^^) as well as 
the renormalization-group equation, ^^^'^^^ shows that Cy appears as a part of the effective 
interactions together with those in the scalar channels, which are responsible for the chi- 
ral symmetry breaking (xSB): The instanton-anti-instanton molecule model gives for the 
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effective interaction between quarks 



2 _ 

l^raolsym = Gmol{l^ + (^'^"^75^)^] 

c 

[(VSt^TmV')' - (VSr"7.75V')'] + ^(V57m75V')'}+>C8, (1-3) 

wliere r" = (r, 1) and Cs denotes tlie color octet part of the interaction, which we shall 
not write down. Near the phase transition, the instanton molecules are polarized in the 
temporal direction, Lorenz invariance is broken, and thus the vector interactions are modified 
as {^^^r-^Y {'^^qF'^Y. Notice that the instanton-induced interaction breaks the Ua(1) 
symmetry. In reality, however, there should also exist UA(l)-symmetric interactions such as 
the one-gluon exchange interaction or its low-energy remnant as 

4l = Gl {{i^Ll.^L? - [i^Lli^Lf] , (1.4) 



where i/jl denotes the left-handed quark field. It is shown using the renormalization group 
equation that the strengths of the UA(l)-symmetric and violating effective interactions are 
of the same order near the Fermi surface. Thus one sees that the vector interaction exists 
together with other chiral invariant terms which are usually used. Therefore, one may say 
that the previous works dealing with the xSB-to-CSC phase transition without incorporating 
the vector interaction Cy are all incomplete, because this interaction may alter the nature 
of the phase transition significantly. 

We shall show in this paper that the inclusion of the vector coupling induces a novel 
interplay between the xSB and CS through the difference of the respective favoring baryon 
densities and changes both the nature of the phase transition and the phase structure in 
the low temperature region drastically *\ The resultant phase diagram and the behavior 
of the chiral and diquark condensates as functions of (T, n) will be found to have a good 
correspondence with those given in two-color QCD on the lattice. It is thus found that our 
simple model gives a possible mechanism underlying the lattice results. 

This paper is organized as follows. In the next section, the Lagrangian to be used 
is introduced. In §3, we shall give the thermodynamic potential and the self-consistency 
condition for the quark condensate and the pairing field. Numerical results are presented 
in §4. The final section is devoted to a summary and concluding remarks. The appendix 
summarizes the effects of the vector interaction on the chiral transition when the CS is not 
incorporated. 

Preliminary results have been reported in Ref. 52) 
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§2. Model 



As a chiral effective model which embodies the vector interaction as well as the usual 
scalar terms driving xSB, we use a simple Nambu-Jona-Lasinio (NJL) model with two flavors 
{Nf — 2) and three colors {Nc — 3), following Ref. 23). The NJL model may be regarded as a 
simplified version of that with instanton-induced interactions and can also be derived using a 
Fierz transformation of the one-gluon exchange interaction with heavy-gluon approximation 
(see 10),ll),53)-55)). This effective model has the merit that it can be used to investigate the 
chiral transition and CS simultaneously, and hence describes their interplay. It was shown^^^ 
that the physical content given with the instanton modeP^^ can be nicely reproduced by the 
simple NJL model with a simple three-momentum cutoff. This means that although there 
are several choices for the high momentum cutoff which mimics the asymptotic freedom, 
the magnitude of the gap is largely determined by the strength of the interaction and is 
insensitive to the form of the momentum cutoff .^'^^ The Lagrangian density thus reads 

C^Co + Cj, (2.1) 

where 

Co = ^{h ■ d - m)jp, (2.2) 
with m being the current quark mass matrix m = diag(m„, m^), and 

jCi^jCs + U + jCc, (2.3) 

with 

Cs = Gs{{iPi'y + {iPt75TiPf}, (2.4) 

and 

Cc = Gc{{iph5r2XAi^^){ip^h5T2XA^) + (^TsA^V^^) (V'^'tsAaV^) } . (2.5) 

jCv is given in (1.1). Here, ip'^"' = Cip'^, with C = 27^7'^ being the charge conjugation 
operator, and T2 and A^'s are the second component of the Pauli matrix representing the 
flavor SU(2)j-, and the antisymmetric Gell-Mann matrices representing the color SU(3)c, 
respectively. The scalar coupling constant G5 = 5.5 GeV~^ and the three momentum cutoff 
A — 631 MeV are chosen so as to reproduce the pion mass = 139 MeV and the pion 
decay constant /^r = 93 MeV with the current quark mass = md — 5.5 MeV;-*^^) we have 
assumed isospin symmetry. It should be noted that the existence of the diquark coupling Gc 
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and the vector coupling Gy do not affect the determination of the pion decay constant and 
the chiral condensate. Although, there are several sources to determine the diquark coupling 
such as the diquark-quark picture of baryons,^^^"^^^'^'''-* the instanton-induced interaction,^^-* 
renormalization group analysis, and so on, we shall take Gc/Gs = 0.6, which accurately 
reproduces the phase diagram obtained with the instanton-induced interaction.^^) As for the 
vector coupling, we vary it as a free parameter in the range of Gy/Gs = — 0.5 to see the 
effect of the vector coupling on the phase diagram. We remark that the vector coupling 
Gy is given by 0.256*5 in the instanton-anti-instanton molecule model^^^ and 0.5Gc in the 
renormalization-group analysis;^'^-''^^-' the range we employ for Gy thus encompasses these 
physical values. 

§3. Thermodynamic Potential and Gap Equations 

In this section, we calculate the thermodynamic potential in the mean-field approximation 
and derive the coupled gap equations for the chiral and diquark condensates. 
The thermodynamic potential f2 is defined by 

12 = -Tin Tr e"^^, (3.1) 

where (5 — 1/T is the inverse temperature and K — H — ^N, with H and N being the 
Hamiltonian and the quark number operator, respectively. The expectation value of the 
quark number is given by 

N,^{N), (3.2) 

where 

{6) = Tr e-^^^-^^O (3.3) 
denotes the statistical average of O. The quark number density is given by 

P, = A^,/V^=(V'7V), (3.4) 

where V denotes the volume of the system, and it is assumed that the vacuum contribution 

to the quark number is subtracted. *•* 

The quark number Ng can be calculated by means of a thermodynamic relation from i? 

as Nq = —df2/dn, and accordingly, pg is obtained from the thermodynamic potential density 

The rotational invariance of the system, which we assume, implies that the spatial component of the 
expectation value {'tpY'^) vanishes. 
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as*) 

P, = (3.5) 

Since quarks have baryon number 1/3, the baryon number density and chemical potential 
are given by = 1/3 • and plb — 3//, respectively, where iso-spin symmetry is assumed. 
We shall use the quark number density Pq and chemical potential pL for the formulation, but 
Pb and plb will be used in the presentation of the numerical results in §4. 

To apply the mean- field approximation (MFA), we first assume that the system has 
a quark-antiquark condensate and a diquark condensate (■0^^75X2 A2'0), where is 

restricted to A2 owing to the color SU(3)c symmetry. In the MFA, K is replaced by 



il^l-i^ • V + (m + Md) - (pt - 2GvPqho\ip + -{A*iP^ij5T2X2ip + h. c.) 



^lG~s + " ^^^^ 



(3.6) 



Here, and A give the dynamically generated quark mass and the gap due to the CS, 
respectively: 

Mij = -2G5(^^), A^-2Gc{^^'tl5T2\2ij)- (3.7) 

We notice here that pt in Kmfa appears in the combination 

pi - 2GvPg = jl. (3.8) 

Thus, the thermodynamic potential ^^mfa in MFA per unit volume is calculated to be 

u{MD,A;T,pi) = Hmfa/V 

_Ml \A\^ 
-AG's^lG^'^''^' 

+sgn{iJ) e_ + e+ + 2riog (l + e'^^^^-)^^-) (l + 6"^^+)} , 

(3.9) 

where Ep = ^Jf^^TM^, = Ep±jl and e± = ^/^^+~\A^, with 



M^m + Mo (3.10) 



This is a familiar procedure in the a-u model.^^^ 
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being the total (constituent) quark mass and sgn{^^) the sign function. Our thermody- 
namic potential reduces to those given in Refs. 23), 56) when Gy = 0. The quark density 
Pq appearing in (3.9) is expressed as a function of the condensates (M£),Z\) through the 
thermodynamical relation (3.5) as 



n(e-)-n(e+) 



tanh 



pe. 



+ ^tanh^ 
e+ 2 



(3.11) 



where n{^±) is the Fermi distribution function: n{^±) — l/(exp{/3^± + 1}). Equation (3.9) 
together with Eq.(3.11) gives the thermodynamic potential u; with the condensates (M^, A) 
being the variational parameters at given {T,ii); their optimal values give the absolute 
minimum of cu. 

The chiral and diquark condensates in the equilibrium state at given {T, /i) should satisfy 
the stationary conditions for the thermodynamic potential, 



duj 



dM, 



D 



dui 



Md 



0, 



(3.12) 



which are reduced to the self-consistency conditions for the two condensates, 

d^p 1 



Md = SGsM 



A = 8GcA 



(27r)3E^ 

d^p 



l-n(e-)-n(eH-) + ^tanh^ 



+ ^tanh^ 
£-4- 2 



— tanh 
e_ 2 



(27r)3 

Here we have utilized the chain rule 



1 , /3e^ 

H tanh 

e+ 2 



D 



dpq 



duj 

A dpq 



(3.13) 
(3.14) 

(3.15) 



dMo A,p, dMo 

and that for the /A-derivative, together with the fact that Eq. (3.11) ensures the relation 

dpq Md,A 



0. 



(3.16) 



In analogy to the BCS theory of the superconductivity, we call Eqs. (3.13) and (3.14) the 
gap equations. Notice, however, that a solution of the gap equations may only give a local 
minimum, or even maximum of the thermodynamic potential, and it is only a candidate of 
the optimal value of the condensates; one must check whether it gives the absolute minimum 
of the thermodynamic potential. 

Prom the structure of the coupled gap equations and the thermodynamic relation (3.11) 
for Pq, one can extract some interesting properties of the condensates (M^?, A) as functions 
of (T, /x) and also of {T,pq). 
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(1) Once the absolute minimum of the thermodynamic potential and, accordingly, the 
optimal value of {Md, ^) are found at given (T, fi), the quark density pq is given by Eq. (3.11). 
The coupled gap equations (3.13) and (3.14) show us that the optimal value of (M/j, Z\) is 
a function of T and jl, and in this way the possible Gy dependence is absorbed into fi. 
Furthermore, if {M£){T, /i), A{T, n)) is a solution of the coupled gap equations with Gy = 0, 
then 

{Mn{T,fl),A{T,fl)) = (Mo(r,/. - 2GvPg), A{T, - 2GvP,)) (3.17) 

is a solution with Gy 7^ 0. Thus, the whole solution as a function of p is shifted toward 
larger p by an amount 2Gypq. 

(2) Next, we shall examine how the solutions of the coupled gap equations behave as 
functions of (T, pg) instead of (T, p) . Let the 0-th order approximation of the condensates 
be given. Then Eq. (3.11) gives /i as a function of (T,pq), i.e., p = p(T,pq). Thus, 
the first-order approximation of the condensates (M£),Z\) is given as the solution to the 
coupled gap equations (3.13) and (3.14), which are only dependent on T and p, not on p 
and Pq, separately; possible Gy dependence is absorbed into p{T,pq). Thus, repeating this 
procedure, one sees that (M/j, A) becomes only a function of T and pq and is independent 
of Gy, because p, through which Gy can affect the formulas, actually only plays the role of 
a dummy variable. Thus we have proved that there is no effect of the vector interaction on 
the behavior of the solution to the coupled gap equation as functions of (T,pq). 

(3) Does it mean that there is no trace of the presence of the vector interaction in the 
phase diagram in the T-p plane? The answer is no. The effect of the vector interaction 
manifests itself in the critical point or line when the transition is first order. In this case, 
there are several solutions to the coupled gap equations (3.13) and (3.14), corresponding to 
the local minima, maxima and even saddle points of lo; notice that these solutions correspond 
to different baryon densities. Since the thermodynamic potential (3.9) is explicitly dependent 
on Gy in a combination with the quark density, the relative magnitudes of the local minima 
change and can be altered with the vector interaction: In Fig. 1, the right (left) figure in 
the upper panel shows the contour map of the thermodynamic potential uj{M£,,A^ with 
Gy/Gs = 0.2 {Gy/Gs = 0) at T = and ps = Pbo = 1035 MeV, which is actually found to 
be the critical point. The thermodynamic potential u}{Md, Z\) as a function of Md at given 
A — Q, 25, 50 and 80 MeV; i.e., the cross sections along the lines shown in the upper panels 
are given in the lower panels, where the solid (dashed) lines denote u) with Gy/Gs = 0.2 
{Gy/Gs = 0.). One clearly sees that the vector interaction increases the thermodynamic 
potential in the small M^, region for every A; notice that the system with smaller Mu is 
at higher density, as discussed in §1. Thus the absolute minimum given with /A ~ 50 MeV 
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Fig. 1. The upper panels show contour maps of the thermodynamic potential uj at {T,^b) = 
(0,1035) MeV with Gy/Gs = (left panel) and Gy/Gs = 0.2 (right panel) in the Md-A 
plane. The difference between the values of u for the adjacent contour lines is 7.5 x 10^ MeV"^. 
The lower four panels are cross sections of u> cut in a plane with A fixed at 0, 25, 50 and 80 
MeV. The solid (dashed) lines represent the Gy/Gg = 0.2 {Gy/Gs = 0) case. 

at small Mr, when Gy — ceases to be even a local minimum with finite Gy/Gs, and the 
local minimum at Mr, ~ 300 MeV with Z\ ~ in turn becomes the unique local, and hence, 
the absolute minimum. Thereby the double-minimum structure disappears, and the first 
order transition is altered to a crossover. In short, the critical temperatures and densities 
at which the transition from one local minimum to the other occurs are strongly affected by 
the vector interaction, and the critical line of the first-order transition in the T-pg plane is 
changed with the vector interaction. 

In passing, we remark that pq cannot be interpreted as a variational parameter with which 
the thermodynamic potential is minimized: Since Eq. (3.11) is obtained by the stationary 
condition Eq. (3.16), one might have imagined that the equilibrium state could be deter- 
mined by searching for the minimum point of the thermodynamic potential with pg being a 
variational parameter together with Md and A^"^^ However, Eq. (3.16) is found to give a 
local maximum of the thermodynamic potential. That is, the absolute minimum of the ther- 
modynamic potential in the M-A-pg space, if it exists, does not give the thermodynamical 
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Fig. 2. (a) The phase diagram in the T-jib plane with Gy = 0. There are four phases; the xSB, 
CSC, Wigner, and coexisting phases. The small panel is an enlargement around the border of 
the xSB and CSC phases at T = 0. The solid line represents the critical line of a first-order 
phase transition, the dashed line a second-order phase transition, and the dot-dashed line a 
crossover, (b) The corresponding phase diagram in the T-pB plane in units of the nuclear 
matter density pQ. There exist mixed phases corresponding to the first-order transitions seen 
in (a). 
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equilibrium state. 

§4. Numerical results and discussion 

In this section, we show the numerical results and discuss the effects of the vector coupling 
on the phase diagram, the T-jib and T-pB dependence of the order parameters. 

4.1. Phase diagram with no vector interaction 

As preliminary to the discussion on the effects of the vector interaction, we first present 
the phase structure without the vector interaction. We shall show that a coexisting phase 
appears where the quarks with dynamically generated mass are color superconducting. This 
is a manifestation of competition between the xSB and CSC phase transition. 

In Fig. 2(a), the phase diagram in the T-ji plane is shown. One can see that there are 
four different phases, i.e. the xSB phase, the normal quark phase, which we call the Wigner 
phase, the CSC phase, and a "coexisting" phase of xSB and CS; as seen from the upper 
small panel, which is an enlargement of the part around the solid line near zero temperature, 
the last phase, in which quarks with dynamically generated mass are color superconducting, 
occupies only a small region in the T-fj, plane near zero temperature with shghtly smaller 
than /iso = 1035 MeV. 

In the figure, the critical line of first- and second-order transitions are represented by 
the solid and dashed lines, respectively; notice that there exists a dashed line in the upper 
small panel. We remark that there are three kinds of first-order transitions: x SB- Wigner, 
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Fig. 3. The order parameters Md and at T = as functions of /LtB with Gy = 0. There are 
discontinuities of the order parameters at fi^o = 1035 MeV. An enlarged figure near the critical 
point is also shown. The gap A is finite even in the region fiB < A*so- 



XSB-CSC and coexisting-CSC transitions. An artificial critical line of the crossover chiral 
transition is also shown by the dash-dotted line on which the dynamical quark mass takes 
the same value as that at the endpoint Md = 186 MeV, so that the crossover critical line 
is connected continuously with the critical line of the first-order transition at the endpoint. 
*^ With this definition of the critical line for the crossover chiral transition, the critical 
temperature at vanishing chemical potential (// = 0) is found to be 187 MeV, which is 
slightly larger than the critical temperature obtained in simulations of lattice QCD with two 
fiavors.^^ 

In accordance with a widely accepted view,^°)'^^) one sees that the chiral transition is 
first-order at low temperatures: The critical line of the first-order transition emerging from 
a point in the zero-temperature line terminates at 

(Te,/XBe) = (47,990) MeV. 

The figure also shows that the phase transition from the CSC to the Wigner phase is second 
order when T is raised in our model, in which the gluon fields are not explicitly included. 
We have found that the xSB phase is transformed into the coexisting phase at very low 
temperatures by a second-order transition when fj, is raised, as shown in the small panel. 

To see more detail of the coexisting phase, we show the /x dependence of Md and A 
at T = in Fig. 3. One can see that Md (A) shows a discontinuous decrease (increase) 
at fj^B = fJ^Bo, which clearly indicates a first-order chiral (CSC) transition at this point. A 
notable point here is that Md has a finite value even in the CSC phase, because of the finite 
current quark mass; notice that Md is proportional to the chiral condensate and not the 
total (constituent) quark mass M — m + Md- Although we do not show the result here, 
we have checked that Md vanishes in the CSC phase in the chiral hmit (nevertheless see 
Fig. 7 for Gy/Gs = 0.2). On the other hand, there is a region in which A becomes finite 

*^ We have followed the criterion used in Ref. 36). 
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in the ^SB phase, which remains the case in the chiral hmit. We have called this phase the 
coexisting phase. 

We notice that a coexisting phase similar to ours was obtained in some previous works. ^'^-''^^^"^^^ 
For example, Rapp et al.^^^ showed that the instanton-anti-instanton molecule model ad- 
mits such a coexisting phase at finite chemical potential, although the phase structure at 
T 7^ was not examined. However, they questioned the robustness of the existence of the co- 
existing phase, because other calculations using a similar NJL-typc chiral model, ^'^^ in which 
the effective scalar coupling constant in our notation is relatively large, did not exhibit 
such a coexisting phase. In fact, we have also checked that if a slightly larger Gs is used, 
the coexisting phase disappears even in our case. We shall show, however, that the vector 
interaction induces a competition between the ^SB and CSC phase transition, and thereby 
the existence of the coexisting phase always becomes possible with a sufficiently large vector 
coupling. **^ 

The phase diagram in the T-pB plane is shown in Fig. 2(b). This phase structure is 
schematically presented in Fig. 4. Corresponding to the three types of first-order transitions 
mentioned above, there exist three mixed phases, which we call 1, 11 and III, respectively: I 
is a mixed phase of the xSB and Wigner phases, while 11 and 111 are mixed phases of the 
xSB and CSC phases, and the coexisting and CSC phases, respectively. We remark that 
various other mixed phases are possible when the CSC phase is incorporated than when it 
is not. 

This ends our investigation of the phase structure without the vector interaction. When 
the vector interaction is included, the phase structure may be changed significantly, which 
we shall show is indeed the case in the next subsection. 

4.2. Phase structure with the vector interaction 

In this subsection, we discuss effects of the vector interaction on the phase structure of 
hot and dense quark matter by varying the vector coupling Gy by hand in the range of 
Gv/Gs = — 0.5. One will see that the vector interaction causes a nontrivial interplay 
between the xSB and CS phase, causing the optimal condensates to greatly fluctuate in a 
combined way. This effect was not elucidated in the previous work.^^^ 

The phase structure in the T-ji plane with GyjGs — 0.2 is shown in Fig. 5(a). The 
phase diagram consists of the xSB, Wigner, CSC and coexisting phases, as in Fig. 2(a). 

*^ In Ref. 56), the full coupled gap equations for Mjj and A were not solved, which is necessary to find 
the coexisting phase. The coexisting phase discussed in Ref. 64) using the NJL model is thermally unstable. 

**^ We remark that if the ratio Gc/Gg in our notation is artificially large, the coexisting phase can be 
also realized in a broad region in the T-ji plane, as shown in the random matrix model^^^ and in the NJL 
model.^3) 
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Fig. 4. A schematic figure accounting for Fig. 2(b). There are three mixed phases in the T-ps 
plane: I is a mixed phase of the xSB and Wigner phases, while II and III are mixed phases of 
the xSB and CSC phases, and the coexisting and CSC phases, respectively. 
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Fig. 5. The phase diagrams with Gy/Gs = 0.2 in the T-jj, plane (a) and T-p plane (b). The solid 
line represents the critical line of a first-order phase transition, the dashed line a second-order 
transition and the dot-dashed line a crossover. 



The dash-dotted line is the contour line at Md — 198 MeV and is supposed to denote the 
critical line of the crossover transition; the solid and dashed lines represent the critical lines of 
the first-order and second-order transitions, respectively, as in Fig. 2(a). The corresponding 
phase diagram in the T-p plane has the three mixed phases I, II and III which are seen in 
Fig. 4, as well as the xSB, Wigner and CSC phases. 
From these figures, the following points are notable: 

(1) The endpoint of the first-order transition moves toward a lower temperature and higher 
chemical potential, 

(Te, iiBe) = (27, 1056) MeV. 

(2) The chiral restoration is moved toward larger p. This is because the gap equations 
(3.13) and (3.14) are functions of T and p, and thus the explicit Gy dependence is absorbed 
into p, as shown in Eq. (3.17). This means that p given by a fixed Md is shifted toward 
larger values as Gy is increased. 

(3) The region of the coexisting phase becomes broader in both T and p directions in the 
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T-/i plane, and hence also in the T-p plane. This feature is determined dominantly by the 
behavior of Md in the xSB phase. As an example, Md together with Z\, as a function of 
/iB at T = is shown in Fig. 6; the same quantities in the chiral limit are shown in Fig. 7. 
One sees that there appears a small region of ^u^, smaller than but near ^bo in which 
{A) shows a gradual decrease (increase); accordingly, finite and A coexist in this region. 
One should notice here that although the coexistence in this sense is realized even when 
A'-B > /^_Bo, as seen in Fig. 6, Md in the chiral limit vanishes identically in this region, while 
the coexistence of Md and A remains for at jib < fJ'Bo as seen in Fig. 7. In fact, this is 
also the case when Gy = 0, as was noted in §4.1. Thus, calling this the "phase coexisting" 
makes sense. Anyway, the gradual change of the order parameters means that the first-order 
transition is weakened. The decrease of also implies that of the total quark mass M, 
leading to a growth of the Fermi surface for a given fig. The larger the Fermi surface, the 
larger the gap A, owing to the BCS mechanism. Thus the region of the coexisting phase 
in the T-fiB plane becomes broader. This feature can be applied to the case for T ^ 0. In 
short, the vector interaction promotes the formation of the coexisting phase. This is one of 
the points which Buballa et al.^^^ did not clarified, because they used a fixed vector coupling. 
*^ It would be interesting to explore the possible correlation between the appearance of the 
coexisting phase and the strength of the effective vector coupling extracted, say, from the 
baryon- number susceptibility,^^^'^^^'^^^ as was done in Ref. 68). 
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Fig. 6. The order parameters Md and A as a function of at T = with Gy/Gs = 0.2. 



The characteristics (1) and (2) of the effects of the vector interaction have been known to 
exist for the chiral transition without the CSC transition incorporated. (An account 
of the phase structure without the CSC transition is presented in Appendix A as a reference.) 

However, when the interplay between the xSB and CS phases enhanced with the vector 

*^ The coexisting phase does not appear in Ref. 48), although a relatively large ratio Gy/Gs = 0.5 in 
our notation is adopted. However, we should also notice that a larger Gs leading to a larger constituent 
quark mass M than ours is used there. This suggests that if the driving force responsible for the xSB phase 
as represented by Gg is strong, a larger ratio Gy /Gs is needed for the realization of the coexisting phase. 
It is worth mentioning in this respect that the coexisting phase is obtained in the two-color QCD on the 
lattice in a robust way.^^ 
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Fig. 7. The order parameters Md and A as functions of at T = in the chiral Umit with 
Gy/Gs = 0.2. The other parameters are shghtly changed so as to reproduce the physical 
quantities in the chiral limit: Gg = 5.01 GeF"^ A = 650 MeV and Gc = S.llGeV"^. One 
can see that the chiral condensate vanishes completely in the CSC phase while it has a finite 
value in the coexisting phase. Thus the phase transition from and to the coexisting phase can 
be unambiguously defined. 



interaction is taken into account, the variation of the phase diagram becomes not so simple 
for larger Gy- In Fig. 8(a), we show the phase diagram in the T-yU plane with Gy/Gs = 0.35. 
It is noteworthy that there appear two endpoints at both sides of the critical line of the first- 
order transition. Accordingly, the coexisting-CSC transition at low temperatures becomes 
a crossover transition. We have checked that the crossover transition becomes second order 
in the chiral limit, and hence a tricritical point appears instead of the endpoint of the first- 
order transition. As far as we know, this is the first time it has been shown that the critical 
line of the first-order transition for the chiral restoration can have another endpoint on the 
low temperature side, implying that the transition from the xSB phase to the CSC phase 
at low temperatures becomes a crossover (second order in the chiral limit). Nevertheless it 
is noteworthy that the two-color QCD on the lattice at nonzero temperature and chemical 
potential gives a similar phase diagram; see Fig. 1 of Ref. 8). Again, the lattice result 
might be interpreted in terms of the effective vector coupling, which deserves exploration 
for the purpose of understanding the underlying physics. The two-endpoint structure does 
not appear in the random matrix model with two colors, in which, however, only the two 
auxiliary fields a ~ {il^il^) and A ~ {ipCij5T2T2ip) are explicitly introduced, but not with the 
vector field. It would be intriguing to study whether such a phase structure can be realized 
in the random matrix model with the incorporation of the vector field as an auxiliary field. 

The corresponding phase diagram in the T-p plane is shown in Fig. 8(b). Its schematic 
phase structure is represented in Fig. 9. The phases II and III correspond to the mixed 
phases of the xSB and CSC phases, and the coexisting and CSC phases, respectively as in 
Fig. 4. Notably, the phase I does not exist anymore. 

To examine the mechanism of the appearance of the two end points in detail, we show 
the thermodynamic potentials in the Mjj-A plane for various T and /i in Fig. 10. In the 
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Fig. 8. The phase diagram with Gy/Gs = 0.35 in the T-/x plane (a) and T-p plane (b). There 
appear two endpoints of the first-order transition. 
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Fig. 9. A schematic figure accounting for Fig. 8(b). The mixed phase I of the xSB and Wigner 
phases does not exist in this case. 



lowest panels, the thermodynamic potential at T = 5 MeV is shown: We see only one local, 
and hence, the absolute minimum point, which varies continuously as fi is increased. This 
implies that the phase transition is a crossover at T = 5 MeV. At higher temperatures, 
however, the thermodynamic potential comes to have two local minima near the critical 
point, as shown in the second and third row panels for T = 12 MeV and T = 15 MeV, 
respectively, and the phase transition becomes first order. At even higher temperatures, the 
double-minimum structure ceases to exist and the thermodynamic potential has only one 
local minimum again, as shown in the uppermost panel for T — 22 MeV, and the phase 
transition again becomes a crossover. 

In our model calculation, the two-end-point structure of the phase diagram appears for 
finite Gy but in a narrow range of Gy/Gs, i-e. 0.33 < Gy/Gs ^ 0.38. We should also 
note that even when the phase transition is first order, the height of the bump between 
the two local minima of the thermodynamic potential per particle is so small that it is 
found to be comparable with or smaller than the temperature. This means that thermal 
fluctuations, which are ignored in the mean-fleld approximation employed in this work, may 
easily destroy the two-end-point structure. What we have found is that the inclusion of the 
vector interaction makes the minimum of the thermodynamic potential shallow in the Mjj-A 
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Fig. 10. The contour of the thermodynamic potential in the Md-^ plane for various values of 
(T, hb) around the critical point of the first-order transition. The difference between the values 
of u for adjacent solid lines is 1.5 x 10^ MeV^. As shown in the second and third row panels, 
there appear two local minima at T = 12 MeV and T = 15 MeV with hb near the critical 

value, which indicates that the phase transition is first order at these temperatures. On the 
other hand, as shown in the bottom and top panels, there always exists only one local, and 
hence, the absolute minimum at T = 5 and 22 MeV, which minimum moves continuously as 
is increased, implying that the phase transition is a crossover. 
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Fig. 11. (a) The phase diagram in the T-/x plane with Gy/Vs = 0.5. (b) The corresponding phase 
diagram in the T-p plane. 



plane, suggesting the significance of fluctuations of the chiral and diquark condensates in a 
combined manner. The incorporation of the thermal fluctuations is beyond the scope of this 
work. 

When we consider a larger value of Gy than 0. 386*5, the first-order transition disappears 
and it is changed completely into a crossover transition. Figure 11(a) shows the phase dia- 
gram in the T-ji plane with Gv/Gg = 0.5 as a typical example in this case. The dashed line 
denotes the second-order transition. The dash-dotted line represents an artificial crossover 
fine on which — 200 MeV. The corresponding phase diagram in the T-p plane is shown 
in Fig. 11(b). One can see that there is no first-order transition, and hence no mixed phase. 
As pointed out in §3, the vector interaction affects the phase diagram in the T-p plane only 
when there is a first-order transition. Therefore, the phase structure in Fig. 11 no longer 
changes after Gy/Gs exceeds 0.38. 

As a nice summary of the effects of the vector interaction on the phase structure of hot 
and/or dense quark matter, we show three-dimensional plots of the dynamical quark mass 
Md and the gap A in the T-fj, and T-p plane in Fig. 12. The thick line represents the critical 
hne of the first-order transition. The dotted points indicate the endpoints. One sees that Md 
decreases more smoothly for larger Gy in the T-p, plane. It is clear that the Gy dependence 
of Md and A in the T-p plane appears only in the critical region of the first-order transition. 

§5. Summary and concluding remarks 

We have investigated effects of the vector coupling on the chiral and color superconducting 
phase transitions at finite density and temperature in a simple Nambu-Jona-Lasinio model 
by focusing on the implication of the density jump accompanied by the chiral transition. We 
have shown that the phase structure is strongly affected by the vector interaction, especially 
near the critical line between the chiral symmetry breaking (xSB) and color superconducting 
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Fig. 12. The left panels show three-dimensional plots of the order parameters Mjj (solid lines) and 
A (dashed lines) as functions of [T^hb) with various values of Gy/Gs, while the right panels 
show them as functions of (T, ps)- The thick line corresponds to the first-order transition and 
the circles denote their endpoints. Notice that the behavior of the order parameters as functions 
of (T,pb) does not depend on Gy, except in the region of the mixed phases, in accordance 
with the discussion given in §3. 
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(CSC) phases: The first-order transition between the ^SB and CSC phases becomes weaker 
as the vector couphng is increased, and there can exist two endpoints of the critical fine of 
the first-order restoration in some range of parameters values; the two endpoints become 
tricritical points in the chiral limit. Our calculation has shown that the repulsive vector in- 
teraction enhances the competition between the xSB and CS phases, leading to a degeneracy 
in the thermodynamic potential in the Mjj-A plane. This implies that there exists gigantic 
fluctuations of the order parameters that appear in a correlated way near the critical re- 
gion, and it suggests the necessity of a theoretical treatment incorporating the fluctuations. 
This is, however, beyond the scope of this work. We have found that the coexisting phase, 
in which the quarks with dynamically generated mass are color-superconducting, appears 
in a wide range of values of jis and T. Here it should be emphasized that it is not yet 
known whether the chiral and the conflnement-deconflnement transitions occur simultane- 
ously at flnite density; hence there may exist quark matter with chiral symmetry breaking. 
We have emphasized that the appearance of such a coexistence phase becomes robust and 
hence universal through the inclusion of the vector interaction. We have also shown that 
the repulsive vector interaction causes the transition from the chirally broken phase to color 
super conducting phase to move toward larger hb- 

Although our analysis is based on a simple model, our flnding that the vector interaction 
enhances the competition between the xSB and CSC phase transitions is universal and should 
be conflrmcd and further studied with more realistic models, including the random matrix 
model and on lattice QCD. In fact, phase structure similar to that found here has been 
obtained in two-color QCD on the lattice, ^^'^^ in which there appear two tricritical points 
related to the chiral and CSC transitions, and also the coexisting phase in a wide range 
of the temperature and chemical potential. It may be possible to intuitively understand 
these results in terms of the effective vector coupling, which can be extracted by calculating 
the baryon- number susceptibility. A random matrix study of the QCD phase diagram^^^'^^^ 
incorporating the vector condensate (i.e., the density) explicitly should be carried out. 

In this work, we have ignored color neutrality in the CSC phase. It is known, 
however, that color neutrality seems to have an only small effect on the onset of CS: When 
the superconducting gap is sufficiently smaller than the Fermi momentum, as is the case in 
a wide region of the T-/i plane in our calculation, the densities of the paired and unpaired 
quarks are close in magnitude. Therefore, the present results obtained for the effects of the 
vector coupling will only slightly change when the color neutrality is taken into account, 
although the two-end-point structure, which is realized through a delicate interplay between 
xSB and CS through the vector coupling, might disappear or persist in the mean-field 
approximation we have employed. 
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Fig. 13. The Gy dependence of the phase diagram for the chiral transition in the T-fx plane. The 
soUd Hne represents the critical line of the first-order transition. The dash-dotted line denotes 
an artificial critical line of the crossover transition, which is determined with the same condition 
as that stated in the text. 

In the two-flavor case which we have treated in this paper, the first two-color states which 
form the color Cooper pairs should have different dynamical mass from the remaining color 
state. Although incorporating such color-dependent dynamical masses^^^'^^^'^^^ is known to 
import only a tiny effect, it should be taken into account for a complete analysis of the effects 
of the vector coupling on the phase boundaries. 

Furthermore, when applying the theory to neutron star phenomena, the charge neutrality 
and the beta equilibrium condition incorporating degenerate neutrinos^^^'^^^ should also be 
taken into account. We have confined our investigation to the two-fiavor case in this work. 
Needless to say, it would be very interesting to examine the effects of the vector interaction 
in the three-fiavor case, and thereby on the color-fiavor locked phase. 

T.Kunihiro thanks David Blaschke for informing him of the work by Buballa et al.^^^ and 
related papers after the completion of this work in May. We are grateful to Michael Buballa 
for pointing out us some misleading statements in the original manuscript with regard to Ref. 
48) and subsequent discussions for elucidating the relation between the present work and 
Ref. 48). M. Kitazawa thanks H. Abuki for communications confirming the precise meaning 
of his master thesis. This work is partially supported by Grants-in-Aid from the Japanese 
Ministry of Education, Science and Culture (No. 12640263 and 14540263). 

Appendix A 

Effects of the Vector Interaction on the Chiral Phase Transition 

In this appendix, we summarize how the chiral transition is affected by the vector interac- 
tion in the case that the CS is not incorporated. Although this problem has been examined 
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Fig. 14. A conceptual diagram accounting for the relation between the total quark mass M = 
m + Mr, and the density pq at given ^u^. The quark density in the equilibrium state becomes 
small for larger M with jiq fixed. 
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Fig. 15. The thermodynamic potential at T = as a function of Md for Gy = 0, 0.1, 0.2, 0.3. For 
smaller Gy, uj has two local minima, reflecting a first-order transition. As Gy is increased, ujd 
at small Mn becomes large and the local minimum at smaller Md disappears. In this way, the 
chiral transition becomes a crossover. 



by some authors, ^^^'^^^ no coherent summary has been given in the literature. 

The phase diagram of the chiral transition in the T-fj, plane is shown in Fig. 13. Here 
we have used the same Lagrangian (2.3) as that used in the text, but with Gq switched off. 
One can see the following features from Fig. 13: 

(i) The chiral restoration is shifted toward larger fj, as Gy is increased. 

(ii) Gy acts to moves the endpoint toward lower T and larger /x. 

(iii) The chiral restoration eventually turns into a crossover transition for large Gy. 

The feature (i) can be understood as follows. The Fermi momentum pp = y^/i^ — 
becomes large (small) for small (large) M, where M = m + Mjj is a constituent (total) quark 
mass, and so does the density p at the fixed /j, (see Fig. 14). Since the vector interaction 
gives rise to a repulsive energy proportional to the density squared, Gyp^, a system with a 
smaller density is favored when Gy is present. Thus one can see when Gy is finite, the larger 
M is favored. We show the thermodynamic potential a; as a function of Mo with various 
Gy in Fig. 15. We see that the thermodynamic potential at small chiral condensate Md 
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increases as Gy increases, owing to the repulsion of the vector interaction. Accordingly, the 
chiral restoration is shifted toward large ^ as Gv is increased. 

Figure 15 also shows that the first-order transition is weakened as Gy is increased: One 
sees from the far left panel that the thermodynamic potential with Gy = has two local 
minima and there exists a bump between these minima. This two-local minima structure 
becomes less prominent and the local minima becomes closer as Gy is increased (see the 
Gy/Gs — 0.2 case (short-dashed line) in the middle panel). Such two- local minima structure 
disappears at Gy — 0.3 for all /xb. 
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